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I investigate the Kazakov-Migdal (KM) model — the Hermitean gauge-invariant 
matrix model on a D-dimensional lattice. I utilize an exact large- iV solution of the 
KM model with a logarithmic potential to examine its critical behavior. I find crit- 
ical lines associated with j str = —1/2 and j s t r = as well as a tri-critical point 
associated with a Kosterlitz-Thouless phase transition. The continuum theories are 
constructed expanding around the critical points. The one associated with ^ str = 
coincides with the standard d = 1 string while the Kosterlitz-Thouless phase transi- 
tion separates it from that with 7 s t r = — 1/2 which is indistinguishable from pure 2D 
gravity for local observables but has a continuum limit for correlators of extended 
Wilson loops at large distances due to a singular behavior of the Itzykson-Zuber 
correlator of the gauge fields. I reexamine the KM model with an arbitrary poten- 
tial in the large-D limit and show that it reduces at large iV to a one-matrix model 
whose potential is determined self-consistently. A relation with discretized random 
surfaces is established via the gauged Potts model which is equivalent to the KM 
model at large iV providing the coordination numbers coincide. 
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1 Introduction 



Matrix models are usually associated with discretized random surfaces or strings in a tri- 
dimensional embedding space and, in particular, with two-dimensional quantum gravity. 
The simplest Hermitean one-matrix model corresponds to pure gravity while a chain of 
Hermitean matrices describes two-dimensional gravity interacting with d<l matter. The 
long-standing problem with this approach is that the stringy phase does not exist above 
the d=l barrier (see Ref. U for a review). 

A natural multi-dimensional extension of this construction is the Kazakov-Migdal 
(KM) model pf which is defined by the partition function 



J n ^n^e^K-^^E^^^^^) (L1) 

{x,y} % 



Here the integration over the gauge field U xy is over the Haar measure on SU(N) at each 
link of a .D-dimensional hypercubic lattice with x labeling its sites and {xy} labeling 
the link from the site x to the neighbor site y. The model (|1 . 1| ) obviously recovers the 
standard open matrix chain if the lattice is just a one-dimensional sequence of points for 
which the gauge field can be absorbed by a unitary transformation of (f) x . 

The KM model was originally introduced in the context of induced lattice gauge 
theories 0. Its remarkable property is an existence of self-consistent scaling solutions 
with nontrivial critical indices [Q for a quartic potential V in ( |1 . 1| ) at large N and any D. 
Some other investigations of a critical behavior of the KM model have been performed 
recently [f|, [5]. There is however a number of problems with the scaling solutions. In 
particular, it is not clear what physical system are they associated with. 

From this point of view it is instructive to look at exact large-iV solutions of the 
KM model and to pass to the continuum by approaching critical points. The only exact 
solutions of the KM model are known for quadratic || and logarithmic @ potentials. The 
spectral density describing the distribution of eigenvalues of a x-independent saddle-point 



matrix, which dominates the path integral ( |1.1| ) at large N, is quite similar in both cases 



to that for a one-matrix model and was obtained by standard methods (see Ref. || for 
a review). The solution of the KM model with the logarithmic potential reduced 0] to 
algebraic equations for the end points of the eigenvalue support — the boundary equations. 

In the present paper I investigate the boundary equations of the KM model with the 
logarithmic potential which explicitly determine the solution! and calculate specific heat 
(or string susceptibility) at large N. The model has a rich critical behavior: a critical 
line associated with j str = —1/2, a critical line associated with j str = and a tri-critical 
point where a Kosterlitz-Thouless phase transition between these two phases occurs. In 
order to construct the continuum theories I perform an expansion around the critical 
points which is quite similar to the one for two-dimensional gravity. The continuum limit 



2 A similar analysis is performed in Ref. 
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associated with 7 str = coincides with the standard d = 1 string (or 2D gravity plus 
critical matter). A Kosterlitz-Thouless phase transition separates this phase from the 
one with ^ str = — 1/2 which is indistinguishable from pure 2D gravity (or a d = string) 
for local observables which live at the same site of the lattice. There is, however, another 
type of observables — extended Wilson loops — for which the continuum limit sets up 
at distances L ~ 1/ \fe (where e characterizes deviation from the critical point) due to a 
singular behavior of the Itzykson-Zuber correlator of the gauge fields. This phenomenon 
occurs in the vicinity of the tri-critical point where the Itzykson-Zuber correlator changes 
its behavior. 

In order to compare with previous results I reexamine the KM model in the large- .D 
limit and show that at large N it reduces for arbitrary V in (|1 . 1| ) and c ~ 1/D to a 
one-matrix model whose potential is determined self-consistently. While this one-matrix 
model has generically 7 str = —1/2, one can obtain 7 str = l/(k + 1) (k > 1) by tuning the 



value of c quite similarly to Refs. [fTO|, |TT|, [I2|j. For the involved logarithmic potential one 



obtains in this way 7 s j r = in agreement with exact solution. 

In order to relate the KM model with discretized random surfaces and strings, I propose 
a matrix model — the gauged Potts model — which is equivalent to the KM model at large 
N providing the coordination numbers coincide. The gauged Potts model has a natural 
connection with discretized random surfaces and is convenient for interpreting the results 
obtained in this paper. The proof of equivalence of the gauged Potts and KM models is 
given via loop equations which reduce at large iV to a one-link equation that is similar to 



the one for the Hermitean two- matrix model |T3| , |14| , |T5|| . This reduction and, therefore, 
the equivalence hold in the strong coupling phase where the vacuum expectation values 
of the closed Wilson loops of the gauge field vanish except for those of vanishing minimal 
area. 

This paper is organized as follows. Sect. 2 is devoted to the description of the KM 
model with the logarithmic potential, its one-cut solution at large N and the calculation 
of string susceptibility. Critical points and the phase structure of the model are obtained. 
In Sect. 3 an expansion around the critical points which specifies the continuum limits 
of the model is performed. The behavior of the Itzykson-Zuber correlator of gauge fields 
in the continuum is studied. In Sect. 4 the large- .D limit of the KM model with an 
arbitrary potential is discussed. In Sect. 5 the gauged Potts model is introduced and its 
equivalence to the KM model at large N is proven. This section also contains a description 
of the large- N solution of these models. In Sect. 6 I discuss the results and some related 
problems for future investigations. Appendix A is devoted to a non-standard behavior 
of the eigenvalue support which leads to •y str = 0. Appendix B contains a proof of the 
convolution formula for the continuum Itzykson-Zuber correlators which is used in Sect. 3. 
Appendix C contains the derivation of the loop equations which are considered in Sect. 5. 
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2 Explicit solution for logarithmic potential 



The distribution of eigenvalues of <p x for the KM model with the logarithmic potential 
coincides at large N with that for a one-matrix model which interpolates between cubic 
and Penner potentials and is solved in this section by the standard technique. The string 
susceptibility of the KM model with the logarithmic potential is calculated at large N. 
It reveals a rich critical behavior: a critical line associated with 7^. = —1/2, a critical 
line associated with ^ str = and a tri-critical point where a Kosterlitz-Thouless phase 
transition occurs. 



2.1 The KM model with logarithmic potential 

Besides the Gaussian case the only exact solution of the KM model (|1.1|) is known for 
the logarithmic potential Jf) i 

V{(j> x ) = -a log (b - <f> x ) - {2D -l)(a + 1) log (a + <f> x ) + [(2D - 1)6 - a]<f> x . (2.1) 



We put for simplicity the coefficient in front of the kinetic term in (|1.1|) c = 1. The 
eigenvalue distribution of the (x- independent) saddle-point configuration coincides for 
the model (|1.1|) with the potential ( |2.1|) and for the Hermitean one-matrix model with 
the potential 

V{4>) = -a log (b - 0) + (a + 1) log (a + 0) - (b + a)(f> , (2.2) 
which is recovered by the potential (|2.1| ) at D = 0. 



The potentials (|2TT| ) and (|2.2|) can be simplified shifting <f> x by a constant value: 

b — a , 
0* = 0* + ~y~ 1 • ( 2 -3) 

The corresponding shift of the kinetic term in the action is 

^ oj '. ni oj = — (f) x u xy (j)yW xy + ~y~ y ^ 0* + ^ <i> y j ■ (2.4) 

Absorbing the last term on the r.h.s. into the new potential, V(4>), and introducing 

P = ^, (2.5) 

one gets 

7(4) = -« log {P ~ k) ~ (2^> - l)(a + 1) log (0 + 4) + 2(£> - 1)^4 (2.6) 

and 

V(k = -a log (p-$) + ( a + 1) log (/3 + 4 - 2 /?0 • (2.7) 



3 The parameter a is related to the original parameters of Ref. M as follows: a = ab + c. 
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Note, that both Eqs. (|2.6|) and ( |2.7| ) can be straightforwardly obtained from Eqs. (|2.1|) 
and ( |2.2j ) substituting a = b = (3. 

The potential ( |2.6| ) depends on two parameters a and f3. However, one more parameter 
7] appears in the perturbative expansion when 

4> x = f)l+(f> ai (2.8) 

and the expansion goes in <p x . Comparing with Eq. (|2.3|) , one identifies r\ with (a — b)/2. 
Therefore, the expression ( |2.1| ) is convenient for the perturbative expansion. 

The perturbative expansion starts from the Gaussian model which is characterized by 
the quadratic potential 

VM = ^ (2.9) 

where 

ml = £ + (2D - 1)- (2.10) 

CI 



as it follows from fl2.1|) in the limit 

a = ab , a~6^oo and <p x ~ 1 (quadratic potential) . (2-11) 

The analogous parameter for the one-matrix model (|2.2|) is 

a b , 

H = t--. (2-12) 
b a 

Solving the quadratic equation ( |2.10| ) for a/b versus m\ and substituting into Eq. ( |2 . 1 2|) , 
one finds 



(D - l)m§ + Dy/mt - 4(2L> - 1) 

(2D^T) (2 ' 13) 



which coincides with the solution of Ref. 0. 

A more narrow region of the parameters than (\2.11 



a = ab , a, b — >• oo , (a — 6) ~ ~ b^ , (2-14) 



is of special interest because the potential (|2.2f) of the associated one-matrix model then 
reduces to a cubic one. 



It is worth mentioning that the normalization of the parameters of the potential Q2.1|) 
is chosen to make (|2.2|) to be D independent. It were be more conventional to have 
instead a D-independent potential V moving the -D-dependence to V. The Gaussian 
formulas (|2.10|) and ( p,13|) show how this can be done. 

For purposes of the perturbative expansion it is useful to restore the hopping parameter 
c in front of the kinetic term, like it enters the exponent in Eq. ( |1 . 1|) . This can be achieved 
by the rescaling 

{(j) x , b} -»• Vc{(f> x , b} , a^—a (2.15) 



at fixed a. This results in the potentials 

V(<j> x ) = -a log (b - <j> x ) - {2D - l)(a + 1) log (a + c<\> x ) + [{2D - l)cb - a]0 x (2.16) 
and 

V{(f)) = - a log {b - 0) + {a + 1) log (a + c0) - (a + c6)0 . (2.17) 
It is easy to see that the potential (|2.17|) of the associated one-matrix model becomes the 



Penner one 16 



Vpenner(<j>) = 01 log (6 - <f)) - a<f) (2.18) 

as c — > 0. 

2.2 One-cut solution 

The one-cut solution for 

where the average is defined with the same measure as in ( II . 1| ) , is given by the general 



formula |T7 

E , = /• ^ H*W(*-*-)(A-* + ) (2 . 20) 
4m (A - u) _ _ x+ ) 

where the ends of the cut, x±, are determined by the asymptotic conditions 

dU V ' {UJ) -0, f duJ . - (2.21) 



c i 2vri {uj — X-){oj — x + ) ^Tci ^ {uj — x^){u) — x + ) 

The contour C\ encircles counterclockwise the cut leaving outside singularities of V'{uj) 
and the pole at A = u so that the integration over uo on the l.h.s. of Eq. ( [2.20|) plays the 
role of a projector picking up negative powers of A. It implies that the branch cut of E\ 
does not pass through singularities (or branch cuts) of V"'(A). This will be the case for 
our meromorphic V^'(A), which has poles at A = —a and A = b, providing they lie outside 
of the cut of E\. 



For V given by ( |2.2| ) the contour integral can easily be calculated taking the residues 
at uj = A, b and —a while the residue at infinity vanishes since falls down as 1/u. One 
gets 



E > 2 



\a\ 1 a + 1 1 



b- A 



yj(b-x-)(b-x+) a + ^ y J(a + x-)(a + x+) i 



(2.22) 



where we use positive numerical values for \J{b — x_){b — x + ) and \J{a + x ){a + x + ) 



which are obtained from the analytic function J (A — X-)(X — x + ) whose cut is depicted 
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—a 



X- 



a) 



-OO 



X. 



X A 



^ b 



b) 




Figure 1: The eigenvalue support of the spectral density (the bold line) and the branch 
cuts of the logarithms (the thin lines): a) for a > 0, b) for a — ► +0, c) for 
— 1 < a < 0, d) for a — > — 1 and e) for a < —1. 
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for positive a in Figs. |T]a), |I]b) and for negative — 1 < a < in Fig. |T]c) . When a — > — 1 
from above, the point x_ approach x + along the real axis as is depicted in Fig. |l]d). For 
a < — 1 the points x_ and x + become complex conjugate as is depicted in Fig. |I|e). 

Some comments concerning Eq. ( |2.22| ) and Fig. [I] are in order. In the Gaussian 
limit ( 2.11 ) Eq. ( |2.22| ) recovers the semicircle distribution of eigenvalues which has the 



support of the type depicted in Fig. [l]a). The analytic function J{X — x)(\ — x+) with 
this branch cut takes at A = b on the value 



^(A-x_)(A-x + )| A=fc = v /(6-x_)(6-x+). (2.23) 

This solution realizes for all a > 0. 

When a — > +0, the end of the cut, x + , approaches b as is depicted in Fig. 0b). For 
a < the boundary equations ( 2.21j) were not have a solution for ( |2.23| ). For this reason 



the branch cut of the square root in Eqs. (|2.20| ) and ( |2.21|) should encircle the point b as 
is depicted in Fig. [T]c) which provides 



- x_)(A - x + )\ x=b = -sj(b - x-)(b - x+) . (2.24) 
Then the boundary equations ( |2.21| ) can be explicitly written as 



lal 1 x_ + x + 
a + + 



(b-x_)(b-x + ) a + b 2 

a + 1 1 x_ + x + 



(2.25) 



y / (a + x_)(a + x+) a + b 2 

They possess a solution both for positive and negative a. 

Eq. ( p.24j ) which is associated with the cut structure of Fig. |l]c) is analogous to that for 
the Penner model |18[] and for the generalized Penner model |19|]. As is shown in the next 
subsection, the boundary equations ( |2.25|) reduce to those for the Penner model under 



the rescaling ( |2.15| ) with c — > at fixed a and the potential ( p. 17] ) smoothly interpolates 
between the Gaussian and Penner ones when a is decreased. 

Eq. ( |2.24j ) always holds when the branch cut of the square root in Eqs. ( |2.20| ) and 
encircles the point b and does not impose further restrictions on the location of the 



cut. It is convenient to take it along the support of the spectral density whose position in 
the complex plane is unambiguously determined by the criterion of Ref . p0 . The proper 
contour for a = —1 is constructed in Appendix A. 

2.3 Solution of the boundary equations 

The boundary equations (p.25 ) can be conveniently rewritten introducing the variable 



a — b 1 x_ + x + , 

z = + t + ^ 2.26 

2 a+b 2 v ' 



8 



and shifting 



- b — a „ b — a „ 6 — a 
{A, rr_, x + | = {A H — , x_ H — , x + H — } 



[2.27) 



2 1 2 2 

which is analogous to (|2.3|) . After this all quantities depend on a and (3 which is defined 
by Eq. fl2.5|) while the dependence on a — b formally disappears from the equations. This 
was already discussed in Subsect. p7T| . In particular, Eq. ( [2.26| ) takes the form 



while Eqs. (|2.25|) read 



at 



(J3-x_)(J3-x + ) 



2(3 



(3 + z 



a + 1 



(j3 + 3;_)(J3 + x+) 



(3-z. 



The product can now be expressed via z by 

^ (a + l) 2 



a 



+ 



2 



_(z + (3f (z-(3)\ 
while z is determined via a and (3 by the 5 th order algebraic equation 

(a + 1) 2 



A/3z-2 



a 



{z-(3f (z + W 
which can be rewritten as 

((3 2 - z 2 ) 2 (4/3z - 2) - (1 + 2a)((3 + zf - Aa 2 (3z = 

Finally, E\ given by Eq. (|2.22 ) takes the simple form 



a 



a + 1 



(3(z + A) 



2 V/3 - A (3 + \J ' v /3 2 -A 2 
Eq. (|2.32| ) is quadratic w.r.t. a. The two solutions read 



[z + (3)((3-z- 



and 



2z' H 2z 2 

1 



It is easy to see from Eqs. ( j2.34| ) and ( |2.35D that 



a + + l = ((3-z)((3 + z- 



2z' 



(3 2 -z 2 







+ 



2z 2 



(2.28) 



;2.29) 



(2.30) 



(2.31) 



;2.32) 



(2.33) 



(2.34) 
(2.35) 



(2.36) 
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so that Eq. Q2.31|) is obviously satisfied. 

Then x± are determined from Eqs. ( |2.28| ) and ( |2.30|) to be 



x± 



l_ ^ - z 2 ) {A(3z - 1) 
2(3 2(3z 



and 



x± 



1 

2~P 



(2.37) 
(2.38) 



for the solutions Q2.34T ) and (|2.35|) , respectively. 

The solution ( [2.38D looks like a meromorphic one and not of the one-cut type. The 
spectral density has, nevertheless, a nontrivial support for this solution at some values 
of the parameters quite similarly for the solution (|2.37|) at a = — l.i The solution (|2.38|) 
can not be obtained, however, from the Gaussian case varying the parameters and the 
susceptibility ( |2.68j ) is strictly divergent at any values of the parameters a and [3. We 
shall consider below for this reason only the solution (|2.37|) denoting a + just as a. 

The equation ( |2.34|) which determines z versus a and f3 is cubic and can be rewritten 
in the standard form 

z z - z([3 2 -a-\) + | = 0. (2.39) 



The discriminant of this equation is non-positive for 

a<(3 2 - 



(2.40) 



When the inequality ( |2.40| ) is satisfied, Eq. ( |2.39p possesses three real solutions. If the 
inequality ( [2.40|) is not satisfied, then Eq. ( |2.39| ) has one real solution. 

As we shall see in a moment the one-cut solution to our one-matrix model exists only 
in the domain of the parameters fl2.40|) where Eq. (|2.39|) has three real solutions. We 
choose the following one of them which is represented in a parametric form as 



z 



_ J (3 2 — a — | cos 
\/3 v 2 



71 

6 



e 

3 







arcsm 



4 V3 



(0> 



4 For the solution ( 2.3£ ) one gets 

a + 1 



E 



A 



o 



/3 - A (3 + X 



■0± 



a 



H - 

/3-A (3 + X 



2(3 



a — 



/3-A /3+A 

-2/3 



(2.41) 



and the support of the spectral density is along the closed contour C in the complex A-plane which is 
determined by the equation 



-2/3A + cdog (/? - A) - (a + 1) log (j3 + A) 



j_ 

'26 







with z given by the solution of the quadratic equation ( |2.35 ). The plus sign should be substituted outside 
of C and the minus sign should be substituted inside C. 
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and < 9 < |. The solution is well-defined in the domain of the parameters a and 
(3 restricted by the inequality ( |2.40| ) which is always satisfied for this solution while the 
equality sign corresponds to sin 9 = 1. This solution is connected to the perturbative 
expansion around the Gaussian model when 9 — > 0. 

In the Gaussian limit (a ~ (3 2 ~ (j3 2 — a) — > oo) one gets from Eqs. ( |2.41| ), ( [2.371 ): 



^-7wh?>' (2 ' 42) 

while (|2.33|) recovers the one for the semi-circle distribution with the parameter 

At = Hj(32- a . (2.43) 



a: 

/i exists for a < (3 2 and vanishes at a = (3 2 which exactly coincides with the criterion 
based on (|2T40| ). 



The limit of the cubic potential is described by the solution (|2.41|) when a —>■ oo in 
the vicinity of f3 2 as (3 — > oo so that 

(3 2 - a ~ /?§ . (2.44) 



No simplifications occurs in this limit in ( |2.41 ) and (2.40) which recover the known results 
for the cubic potential f2lfl. 



In the Penner limit (/3 —>■ oo, a ~ 1) Eqs. ( ^.41| ) and ( |2.37| ) yields 



>=f>-^r, !^i± = v^i (2 45) 

and Eq. (|2.33| ) recovers the solution for the Penner model. The inequality (|2.40|) is always 
satisfied. 

The proper regions of the parameters are depicted in Fig. The bold line which 
starts at a — (3 — corresponds to the equality sign in ( |2.40|) . The solution (|2.41| ) is 
well-defined below this line. 



2.4 The critical behavior 

The critical behavior of our model can emerge when: 

i) The spectral density ceases to be positive at the interval [ for the one- 
matrix models with a polynomial potential. 

ii) Either x_ approaches — (3 or x + approaches (3. 



iii) X- approaches x + like for the Penner model [^2 . 
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Figure 2: The phase diagram of the Kazakov-Migdal model with the potential (2.1). 

The one-cut solution realizes for a < a c . The critical lines a = a c and 
a = — 1 correspond to 7 s j r = —1/2 and 7 str = 0, respectively, while the tri- 
critical point j3 = 1/2, a = — 1 is associated with the Kosterlitz-Thouless 
phase transition. 
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Since we start from the Gaussian solution, we shall verify that a phase transition of 
the type i), which restricts the one-cut solution, does not occur before those of the type 
ii) or in). The proper criterion can be obtained from Eq. ( |2.33| ) where the factor z + A 



should vanish for some A G for this phase transition to occur. The equation for 

the critical points associated with x^ = —z can be easily obtained from Eq. ( [2. 37 ). It has 
the following solution for the critical value of a at given (3 




x_ = -z for a c = (3 2 - 3 y - \ (2.46) 



which exactly coincides with equality sign in fl2.40|) when the discriminant of the cu- 
bic equation ( |2.39|) vanishes. For this reason the inequality X- > —z which means the 



existence of the one-cut solution (|2.33|) is equivalent to (|2.40| ). This is why the expres- 



sion ( p.41|) fully describes our one-cut solution of the one-matrix model. 

It is easy to understand why the critical value (|2.46|) is associated with the vanishing 
of the discriminant. We can identically rewrite Eq. (|2.39|) as 




a c -a= z- T U + - T ( 2 -47) 



which looks like the genus zero string equation in polynomial one-matrix models ]23[. 
For the critical behavior to occur, two roots of the cubic equation should coincide which 
happens when the discriminant vanishes. The critical value of z extracted from Eq. ( |2.47| ) 
is 



(2.48) 



and near the critical point 



a c - a = (z - z c f (l + 2 -^\ w 3(z - z c f (2.49) 

in the full analogy to one-matrix models with polynomial potentials. 

It follows also from the above formulas that z given by Eq. (|2.41 ) is a monotone 



function of both a and j3 in the region where the inequality (|2.40| ) is satisfied since 

da 2 . q o. . 

Tz = ^ {zl ' z ) (2 ' 50) 

and 

d A = z * ~ Z Q (2 51) 



while z > z c in the domain ( 2.40Q . The expression under the square root in the latter 



formula is always non- negative for a > — 1. 
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Another set of singular points is when x + = (3. From Eq. (|2.37|) we find that this 
happens for z given by Eq. ( |2.41| ) only when a = 0: 



x + = (3 for a = , (2.52) 
while z is expressed via (3 along the line a = by 



z = ^ + ^-\ fOT « = °- ( 2 - 53 ) 

Eq. (|2~53|) is applicable for (3 > (3 C = y/2. The value (3 C = V% satisfies Eq. at a = 0. 

Therefore, the critical line ( |2.52| ) which is depicted in Fig. [| by the solid line terminates 
at the intersection with the critical line fl2.46| ). The tri-critical point is given by 

a c = 0, (3 C = V2. (2.54) 



Analogously, x_ and x + coincide for the solution (|2.41|) when a = — 1: 



x_ = x + for a = —1 , (2.55) 

while z is expressed via f3 along the line a = — 1 by 

z = p for a = -l. (2.56) 

Eq. (P3BD is applicable for f3 > (3 C = 1/2. The value /3 C = 1/2 satisfies Eq. at 
a = — 1. Therefore, the critical line (|2.55|) which is depicted in Fig. ||] by the solid line 
terminates at the intersection with the critical line ( j2.46j ). The corresponding tri-critical 
point reads 

a e = -l, (3 C = 1 -. (2.57) 
Note that the critical lines ( j2.52| ) and ( |2.55| ) lie in the allowed region below the line ( |2.46| ) . 



2.5 Calculation of susceptibility 

The critical behavior of the model ( |1.1| ) with the potential ( p.l| ) is characterized by the 
susceptibility 

where Vol. stands for the volume of the system 



Vol. = J2 1 • ( 2 - 59 ) 



By differentiating ( |1.1| ) one gets explicitly 



x = ~L ( I [iog {b ~~ 0x) + {2D ~ i} iog (a + ^ 



fc 2^ [1 ° g {b ~ U) + (2j ° " 1} bg (a + U) ^ (2 ' 60) 



iCi 2iri 
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where 



E,„ = 



da 



(2.61) 



Compressing the contour C\ in Eq. (|2.60|) to the cuts of the logarithms, one gets 

rb 



X 



dtE t + {2D-l) 



dtE t 



This formula can alternatively be derived from (|2.60|) using the identities 



d I tr 



da \ N 



log (b - (f) x ) 



d 

da 



dt 



tr 



and 



d I tr \ d 

-(-log(a + 0,)) = - 



- Q , / tr 
dt ( — 



1 ^ 



N (t — ^) 



dt E t 



dt E t 



(2.62) 



(2.63) 



;2.64) 



To calculate \ in genus zero, we use the formula 
E\ 



dtu V'(uu) - x-)(u - x+) 
Ci 4ttz (A - w) v /(A-x_)(A-a; + ) 



(2.65) 



which can be proven by a direct differentiation of Eqs. (|2.20|) and (|2.21|) and holds for 
the derivative w.r.t. any parameter of the potential V, in particular w.r.t. a. For our 
logarithmic potential ( |2.2| ) one gets explicitly 



b — A a + A 



+ 



yJ(X-X-)(X-X+) 



x 



{b — X-)(b — x + ) J {a + x_)(a + x + ) 



X-b 



+ 



A + a 



(2.66) 



Here and below the plus sign in front of the first term in the square brackets is associated 
with a > while the minus sign corresponds to a < 0. The result of Ref. [B3 for the 
generalized Penner model is recovered by this formula as a — > oo. 

The integral in Eq. ( J2.62p with E t given by fl2.66| ) is easily calculable using the following 
formula for an indefinite integral 

, 2 ^ 



dtE f = - Ioe 



(y/ (t - x-)(a + x + ) + yj(t — x + )(a + x~-j) 

(yj(t - X-)(b - X+) ± yf(t - X + ){b - X-)Y 



(2.67) 



which yields for the susceptibility in genus zero 
Xo = (£>-!) log 



(a + X-){b — x + ) 



± 



(a + x + )(b — X-) 



D\og 



(a + x-) 



\ (a + x + ) \ (a + 



4 l\ (a + x + ){b — xJ) \ (a + xJ){b — x + ) 



(a + x + ) 




\(.t- 



(6-ar_) 



■. (2.68) 
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As before the positive sign in ± corresponds to a > while the minus sign should be 
substituted for a < 0. As is already mentioned in the previous subsection, the a < case 
can be obtained from the a > case by changing the sign of yjb — xJ){b — x+). 

Let us briefly discuss some properties of the expression ( |2.68| ) which will be used in 
the next section for studying the continuum limits of our model. If a > the contour 
of integrations over t, which coincide with the branch cuts of the logarithms depicted in 
Fig. |l|a), can be moved along the complex plane. In particular, one can integrate in the 
first term on the r.h.s. of Eq. (|2.62|) from b to +oo along the positive real axis. The result 
is the same and is not singular at x_ = x + if one substitutes the plus sign in Eq. ( p.68| ). 

On the contrary, the points x_ and pinch the integration contour for a < as 

is depicted in Fig. |l|d) for a = —1. This results in a logarithmic singularity of ( |2.6S| ) at 
£_ = x + when one substitutes the minus sign in Eq. Q2.68Q : 

Xo ~ — log (x- — x + ) 2 . (2.69) 

An analogous logarithmic singularity emerges when x + ~ b as is depicted in Fig. [T]b) at 



a 



0: 



Note, that 



in the Penner limit Q2.15|) with 



Xo « log (b - x+) . 

Xo = -log(l + a) + loga 
-> when x± are given by 

, a + 2 , 2 / - 

x± = b ± — v ol + 1 . 

a a 



(2.70) 
(2.71) 

(2.72) 



Eq. (p. 71 ) coincides with the susceptibility for the Penner model [22] and recovers the 
singularities Q2.69|) and (|2.70|) . 

Having the explicit formula ( P-68| ) for xo-, we can find ° u t which 7 str is associated with 
each type of the critical behavior ( [2.46] ), ( 2.52 ) and (|2.55| ). Along the line (|2.46 ) where 
Xo is not singular and equals to some value Xo j one S e ^ s 

Xo-xZ~{x-~x c -). (2.73) 

Since (x_ — x c _) ~ (a c — a) 1 / 2 near the critical line ( |2.46 ), one obtains 7 str = —1/2. 

Near the critical lines (|2.52j ) and ( p.55| ), where xo is given by Eqs. Q2.69Q and (12.701) , 
one gets 7 str = 0. While (|2.69|) is positive, ( |2.70D is negative. For this reason we shall 
not consider the continuum limit associated with the critical line (|2.52 ) as well as the 



tri-critical point (|2.54j) since it corresponds to a negative susceptibility. The susceptibility 
in the vicinity of the tri-critical point Q2.57|) is considered in Subsect. 3.3. 



3 The continuum limits 



The continuum theories are obtained by an expansion near the critical points. The contin- 
uum limit associated with 7^. = describes in a standard way 2D gravity plus ID critical 
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matter. The continuum limit associated with ^ str = —1/2 looks like pure 2D gravity for 
local observables which are defined at the same site of the lattice. For another type of ob- 
servables — extended Wilson loops — the continuum limit sets up at distances L ~ 1/ ^fe 
with e being a deviation from the critical point. This is due to a singular behavior of 
the Itzykson-Zuber correlator of gauge fields for this phase. A Kosterlitz-Thouless phase 
transition separates the two phases. 



3.1 lsrt = -1/2 

In order to find out what kind of continuum theory is associated with the critical behav- 
ior ( |2.46|) , let us expand all quantities near the edge singularity of the spectral density 



substituting 

z = z c + eVX, Zc= (jj t 3 - 1 ) 

where e — > and A is to be identified with the cosmological constant of 2D gravity. From 
Eqs. (ET47D , one gets 



a c — a = 3e 2 A , x_ = —z c + 2s\fX , x + = 3z c — — . (3.2) 

Introducing the continuum momentum variable, £, by 

A = —z c + e£ (3.3) 

one gets from Eq. (|2.33|) 

E x = ~ + e£- e 3/2 U + VX) ^i+2~7K 2 + 0{e 2 ) . (3.4) 



The last term on the r.h.s. determines the continuum spectral density 

p c (0 = - k + VA) ^-2VX (3.5) 



7T 



and, therefore, all continuum correlators of the trace of powers of the (renormalized) field 
$(x) at some point x, which is the standard set of observables of 2D gravity. 

While the gravitational part of the system is continuous, this does not necessarily 
mean that matter becomes critical. An example is the Ising model on a random lattice 
where it is easy to construct the 7^ = — 1/2 behavior which is associated with continuum 
2D gravity and non-critical matter while matter becomes critical at a tri-critical point 
changing the value of the string susceptibility to 7 sir = —1/3. 

A direct way to verify whether matter becomes critical at a given fixed point is to 
investigate observables which are associated with extended objects — the open-loop av- 
erages 
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where C xy goes from x to y along some path on a D-dimensional lattice and the average 
is w.r.t. the same measure as in flLip. G v \(C xy ) is symmetric in v and A due to invariance 
of the Haar measure, dll, under the transformation U — > IP. 

The averages ( |3.6| ) depend |2|| at large only on the algebraic length L(C xy ) of 
the contour {i.e. the one after contracting backtrackings) and can be calculated [E?5H 
providing C (z/, A) — the one-link Itzykson-Zuber correlator of the gauge fields — is known. 
The latter is expressed via the double discontinuity of 

G vX = G uX (l) (3.7) 

across the cut: 

where the continuous and discontinuous in v parts of G uX are defined by 

G( u +io)\ — G( 



Disc,, G 



v X 



2i 



Cont, G v \ = G(u+l0)X ± G( -° )A (3.9) 

so that for a real A outside of the cut (cuts) Disc„ G u \ coincides with the imaginary part 
and ContyG^A coincides with the real part. In particular, Disc^ E v = —mp(y). 

For the potential Q2.1| ) G u \ reads |7j 

_ (b-u)E v -(a + \)E x + l 
UvX ~ (\ + a + E„)(b-v)-a • [6 Uj 

The few lower terms of the expansion of G v \ in e are 

GvX = iifTT " — —\ — MTT) + ° {£) (3 ' n) 



which is similar to the one in a two- matrix model |[L5| . The double discontinuity of the 
first two terms on the r.h.s. vanishes so that one needs the term at least 0(e). 

For this reason we start directly from the exact C(v, A) which for the potential (|2.1|) 
reads [7j 

C ^ A ) = T h ux ~t4t7T r~y\ (3-12) 

(6-i/)[A-r+(i/)][A-r_(i/)] 

where 

The r.h.s. of Eq. ( |3.12 ) is symmetric in v and A due to the relation 

r±(r T (A)) = A (3.14) 

which is satisfied by the solution (|2.22|) [|J. Analogously to Ref. [[§], it is easy to see that 
Eq. ( |3.14| ) is satisfied by our explicit continuum formulas. 
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To calculate the continuum limit of ( |3.12| ), we expand 

r±(A) = -z c + ei- z jff_ 1} + 0(e 3 ) ± mp(A) . (3.15) 
We put here A = and keep terms of order e 2 which we shall need below. One gets 



1 (4g| -!)(! + e^j^j l+s—^ 
e> {^1 + 1) & - ix y + 



c W A ) - 3773 i in 77 . 2 e ^g A «,+eir~ • ^' lb) 



This expression comes from the 0(e) term in the expansion ( [3 -11 ) since the denominator 
in dp) is 0{e 3 ). 

In the continuum limit e = 0we get from Eq. ( 3.16Q 

C c (x,y)oc- — — (3.17) 
which coincides with C(x,y) for the quadratic potential P5|i 

C(x,y) = - - ^^_ 4 ^ ^ y 2 ^7j (quadratic potential) (3.18) 

at // = 0. While for the quadratic potential p = is possible only at £) = 1 where it is 
associated with the naive continuum limit |J , Eq. ( |3.17[ ) holds in our case at any D along 
the line ( gjgg) . 

The expression (|3.18| ) for the Gaussian C(x,y) possesses a remarkable convolution 
property |26| — keeps its functional structure when one combines the path of the length 
L\ + L2 from the paths of the lengths L\ and L 2 . Then 

C(x, y; U + L 2 ) = J dt p(t) C(x,t; LjC^y; L 2 ) (3.19) 

where 

C{y, A; L) = Disc, Disc A G vX {C xy ) . (3.20) 

7i z p(u)p(A) 

Observables which are associated with matter reveal singularities at p — and become 
finite after an appropriate renormalization. 

Since our expression ( [3.17 ) looks like the Gaussian one for p = 0, one might expect a 
similar property of matter correlators. However, when Q3.17 ) is substituted into Eq. ( |3.19|) , 
the integral is divergent at t — x or t — y. To regularize, we keep 0(e) term in the 
denominator of ( |3.16| ). Now the integral is 0(e~ 1 ^ 2 ) which exactly cancels an extra 0(e 1 ^ 2 ) 
factor which emerges since each of two C in the integral Q3.19| ) is proportional to 1/e 2 
(according to Eq. ( |3.16| )) and the measure is proportional to e 5 ^ 2 . Therefore, one gets 



5 An analogous formula for the two-matrix model has been earlier obtained in Ref. | H | . 
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and the functional structure of C c is preserved. The coefficient of proportionality is 
calculated below. 

An interesting question is whether Eq. ( |3.19| ) with p given by the continuum for- 
mula ( |3.5| ) can have a solution with the same functional structure for macroscopic loops 
as well. A solution to this equation for A = is given by 



C c (x,y; y/u) 



2y/Tl 



(x — y) 2 + 2u(x + y)xy + u 2 x 2 y 2 
which obeys the following convolution property 

1 f°° 3 

dt t*C c (x, t; y/u)C e (t, y; y/v) = C c (x, y,\/u+ y/v) . 



7T JO 



(3.22) 



(3.23) 



This formula is proven in Appendix B. 

It is easy to see that the solution (|3.22j ) satisfies the initial condition (|3.17 ) for u — > 0. 
To find the relation between u and the length L, let us first rescale 



u — > u 



e 4 -l) 



(3.24) 



2 ' 2 

Now the denominator in ( 3.22 ) coincides to order e with the one in ( |3.16|) providing 

u = £ - . (3.25) 

It follows then from Eq. ( 3.23Q that 



u = L — 



(3.26) 



1 

L ~ — -= 

>e 



where the length L is measured in the lattice units. 

One sees from Eqs. fl3.22j ) and (|3.26| ) that the continuum limit of extended correlators 
is reached at distances i 

(3.27) 

rather than ~ 1/e as it might be naively expected. Therefore, a nontrivial scale dimension 
of the matter field is developed at macroscopic distances. 

The knowledge of the L-dependence of C c allows us to calculate correlators of extended 
objects. The simplest one is that of the adjoint Wilson loop 



W A (C) = ^ 
Substituting ( p.22|) , one gets 



£ 17(C) 



1 



i Jdtp(t)C(t,t;L). (3.28) 



W A {C) 



dt t 3 ' 2 



2y/u 

E7T Jo 4ut 3 + U 2 t 4 



1 

'Te 



(3.29) 
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Since the integral is divergent as t — > 0, the result is obtained by an analytic continuation. 
This divergence might be related to the fact that we put A = 0. 

Note, that the correlator ( |3.29| ) is divergent when e — > similarly to the one for the 
Gaussian case [^6j as \i — > 0. One concludes, therefore, that matter could become critical 
at distances ( |3.27| ). This mechanism does not work for the two- matrix model which is 
associated with D = 1/2 in the above formulas. In order to have L — > oo as e — > 0, one 
needs D > 1. For D = 1 one should choose a closed matrix chain since otherwise the 
gauge field can be absorbed by a gauge transformation of <j) x . 

Analogously to ( |3.29 ) one can calculate the more general correlator 



A-0 ; 



-U{C X 



p(t)C(t,t;L) 



1 



A-i 



In the continuum limit one gets 



e 2 7r 



dt t 3 / 2 



(Aut 3 + u 2 t 4 ) (e - t) 




_ -3/ 2 _ ^ 



3/2 



1 + 



This expression is also divergent when e —>■ 0. 

Finally, the continuum part of G U \(C) is given by 



1 /- 00 

Gtf(C) = 2y/eu— / dxdy- 

IX J ( 



1 



(x — y) 2 + 2u(x + y)xy + u 2 x 2 y 2 (£ — x)(^ — ?/) 

(-C) 3/2 (-0 3/2 



2V£M 



(c - o 2 + mc + oce + u 2 ee 



+ 



(3.30) 



(3.31) 



(3.32) 



This expression is indeed 0(e) for L ~ 1 in agreement with ( |3.11| ) while it is 0{\/e) for 
L ~ lUfe. 



3.2 7 s t r = 



Another continuum limit is associated with the critical line a = — 1. One substitutes 
the expansion near the edge singularity which is quite similar to the one for the Penner 
model [18]. Introducing the cosmological constant, A, by 



a: 



-l + e 2 A 4(3 2 ( 3.33) 
where an extra factor is inserted for a latter convenience, we get from Eq. ( J2.41 ) 

z = P- (3.34) 

and from Eq. ( |2.37| ) 



x± = P-^±2eVA. 



(3.35) 
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Defining the continuum momentum £ by 



one gets from Eq. (|2.33|) 



A 



-2(3 - 2(3 2 e£ - eJi 2 - 4A 2(3 2 + 0(e< 



which determines the continuum spectral density to be 

1 



Pc(0 = -V£ 2 ~ 4A - 

7T V 

The one-link Itzykson-Zuber correlator ( 3.12] ) reads 

A(3 2 



C(u,X) 



+ 0(e) . 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



4/3 2 - 1 

This expression does not depend on the continuum momenta £„ and £a similar to the 
Gaussian expression ( 3.18 ) for x, y <C 1. Thus, nothing special happens with the Itzykson- 
Zuber correlators in the 7 sir = continuum limit far away from the tri-critical point in 
contrast to the 7 sir = —1/2 case. 



3.3 Kosterlitz— Thouless phase transition 



The above formulas are not applicable in the vicinity of the tri-critical point (3 C = 1/2, 
a c = — 1 where the continuum system undergoes a Kosterlitz-Thouless phase transition, 



which was previously studied for a closed matrix chain |27| , |28fl , between the phases with 
Istr = — 1/2 and 7 sir = 0. This domain is most interesting and should be treated sepa- 
rately. 



Let us expand 



P = g + ^ 



a 



■1 + 5a 



1 . 
- + Sz 



A 



Then for the solution ( 2.41]) one gets 



5z = -5(3+ J -(5{3) 2 - -5a 



9 



while the inequality ( |2.40|) reduces to 



8a < -(5(3) 2 



(3.40) 



(3.41) 



(3.42) 



It is convenient to parametrize the vicinity of the tri-critical point by the lines 

5a = (3k + 1)(1 - k)(5(3) 2 (3.43) 
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which obey the inequality (|3.42|) for any k. We choose 1/3 < k < 1 where the factor in 
Eq. ( |3.43| ) is a monotone function of k. The cases k = 1/3 and k = 1 are associated with 
the critical lines ( |2.46| ) and (|2.55|) , respectively. Then, ( |3.41| ) is rewritten as 



8z = n8(3 



(3.44) 



while Eq. (|2.37|) gives 



x± 



- + (2 + k)8(3 ± 2^2(1 - k 2 )8(3 . 



(3.45) 



Eqs. ( |3~43D , (^4§, (|3~43p recover Eqs. 



3) or Eqs. (|3~33|) , (CT) (|3735|) if 



1 ey/A 
3 + ^5/T 



or 



1 - 



e 2 A 



(3.46) 
(3.47) 



4(5/3) 2 

with e -C 5/3, respectively. 

The susceptibility ( |2.68|) can easily be expressed via 5(3 and k near the tri-critical 
point: 



Xo 



log 



K 



1 + 3k 



-2D \og(V25{3). 



(3.48) 



The first term on the r.h.s. recovers Eq. ( |2.69| ) in the limit fl3.47| ) while the second one 
which becomes singular when 8(3 — > is a new type of singularity which appears only at 
the tri-critical point. Notice that this second term is D dependent. Using Eq. Q3.43D one 
can rewrite ( [3.48D alternatively via 8a. 

The continuum spectral density can be obtained from ( [2.33| ) by substituting the ex- 
pansions (|3.44j) and ( |3.45|) which gives 



Pc(8X) = ^(8\- (2 + k)8(3)>- 8(1 -K*)(8(3y 



k8(3 + 8X 
8(3 + 8X 



(3.49) 



For k > 1/3 the zeros both of the numerator and of the denominator lie outside of the 
eigenvalue support. In the limits ( |3.46|) and (|3.47| ) this expression recover Eqs. (|3.5| ) and 
( p.38| ), respectively. 

The one-link Itzykzon-Zuber correlator (|3.12|) near the tri-critical point reads 



C(Su, 8Xj) 



{8\ + 8(3){8v + 8{3) 



(3.50) 



where 



V(8u, 8X) = 8\8v(8\ + 8u) + ((8\) 2 + {8uf)8(3 
+k{2 - 3k)(5A + 8u){8(3) 2 + 2k 2 {1 - 2k){8(3) 3 . 



(3.51) 
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This expression is ~ 1/5/3 in agreement with ( [3.16|) and (|3.39|) which are recovered when 
one substitutes 

5X- 



and 



respectively. 



S\ = 35(3 + e£ . 



4 Critical scaling in the large- D limit 



(3.52) 
(3.53) 



The Itzykson-Zuber integral 



dU e Nctr <f>* u <t>y u \ 



(4.1) 



which enters the partition functions ( |1 . 1|) , can easily be calculated at small c (i.e. at strong 
coupling) : 



c 2 N 2 



log iy> x , <f) y ] = ctr^tr <f) y -t 



tr 2 / tr 



N 



N 



tr 2 / tr 



N 



N 



+0(c 3 ). (4.2) 



Let us consider the large- D limit of the partition function. Assuming that V(<j> x ) ~ 1 
as D — ► oo,i we see that the kinetic term is of order D (i.e. of the same order as the 
potential) if 

(4.3) 



1 

D 



The Itzykson-Zuber integral coincides in this limit with the first term of the expan- 
sion ( |4.2[ ) since the higher terms are suppressed. Therefore, we can write down the 
partition function (|1.1|) as 



II' 



(4.4) 



Further simplification occurs in the large- N limit when we can replace one trace in the 
product of two traces in the exponent in ( |4.4| ) by the average value due to factorization. 
One arrives, hence, to the one-matrix model whose potential V(<f>) is determined self- 
consistently from the equation 



V-(0)=V(0)-2cI>(!^ 0. 



(4.5) 



3 We differ at this point from Ref. B where the large-D limit was considered for c ~ 1 and V((j> x ) ~ D 



so that the solution with the minus sign in front of the square root in Eq. (2.13) was chosen in order for 
fj, to be ~ 1. Contrary to the statement of Ref. H] that there is no real scaling solutions at large D, we 
do have them in the limit ((0]) perturbing the Gaussian solution ( 2.1 3| ) with the plus sign which agrees 
with the large-mass expansion. 
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We have assumed that the potential V(4>) is not symmetric so that (tr0) 7^ 0. If V(<p) is 
symmetric {V(<p) = V(—4>)) and (tr0) = 0, then one should keep the second term on the 
r.h.s. of Eq. ( |4.2| ) which yields 



V((j)) = 1/(0) - c 2 D ( ^ 2 \ 2 (symmetric potential) (4.6) 
\ N / v 

and 1 

(symmetric potential) . (4-7) 



Eq. ( |4.6| ) is analogous to that for the one-matrix model with a symmetric potential 
which involves (tr0 2 ) 2 |10|| . An equation of the type ( |4.5| ) appears [II, [T2|] for a non- 
symmetric potential. Quite similar to these papers we shall see in a moment that the KM 
model at large D admits scaling solutions with 7 sir > 0. 

Let us analyze to this aim Eq. (^4.5| ). For the one-cut solution of the Hermitean one- 
matrix model with the potential V one gets 

tr \ r dtu V\uj)uj 2 x_ + x + 



N lv Jci ^iJ(u-x-)(u-x+) 2 



ci ^V'(u) y/(u-x.)(u-x + ) + . (4.8) 



Eq. ( |4.5|) can then be written as the following equation for V: 

duj r — X- + x + 



V(X) = V(X) - 2cD 



,V'(uj) ^J(uj- x-)(u - x+) + 



A. (4.9) 



Differentiating w.r.t. a parameter of the potential which is associated with the cosmolog- 
ical constant and w.r.t. A, we get finally 

dui a 



V\X) = V\\)-2cD j-V'(lu) y/(uj-x-)(uj-x + ). (4.10) 

Let us identify the cosmological constant with g\ — the coupling in front of the linear 
term of the potential. For the susceptibility in genus zero one gets 

*,=/ :^vmk=- 9i { -^^ (4.11) 

JCx 2m 16 

while Eq. (|4~T0| ) yields 

~9i = - n L • ( 4 - 12 ) 

1 + cD (x ~ * +) 

To obtain the critical behavior, we expand near x_ = x c _ which gives for a k-th 
multicritical point of the one-matrix model [p^[ : 

x_ - x c _ ~ {g{ - ~ gi ) 1/k ■ (4-13) 
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Under normal circumstances when Eq. ( |2.73| ) holds, one gets from ( |4.13| ) 

Xo 

so that 7 str = — 1/k since 



This is noi the case, however, for 



when the denominator in Eq. ( |4.12j) vanishes. At this point one has 

9i ~ 

so that 



(rf-<&)~(#-ft) (Jfc+1)/fc ~(*-- 
For the susceptibility (|2.60|) one gets at the point ( |2.46|) 

xiy^igl-gi) 



\k+l 



-l/(fc+l) 



(4.14) 
(4.15) 

(4.16) 

(4.17) 
(4.18) 

(4.19) 



Xo ~ [X- -x^) 
which is associated with 7 sir = l/(k+ 1). 

The formula ( |4.5D , which describes the reduction of the KM model to an one-matrix 
model at large N in the large- .D limit, can be explicitly verified for the potential ( |2.16| ) 
when the exact solution is known at any D. As D —>■ oo with c —>■ according to Eq. ( fl.3| ). 
the potential V of the one-matrix model is given by the Penner potential ( |2.18| ). An 
explicit calculation which is based on Eq. ( 2.22j ) yields 

\N7v a 
while the large-D limit of the potential ( p. 16 ) is 



(4.20) 



V{<f>) -> -a log (b - 0) - a + 2cD 



(a + 1) 



(4.21) 



where (|4.3| ) is used. It is easy to see now that Eq. ( f4.5|) is satisfied. 

While the reduction to a one-matrix model with the Penner potential holds for the KM 
model with the potential ( 2.16|) in the large- .D limit, the only possible scaling behavior 
is with 'jstr = in a perfect agreement with the results of Sect. 2. This seems to be a 
limiting case of j str = l/(k + l) which appear from the critical behavior with j str = —1/k 
of the one-matrix model with the potential V. 

Let us note finally that nonvanishing results for continuum correlators can be obtained 
in the large- .D limit only for these of operators living at the same lattice site while the 
Itzykson-Zuber correlator for a contour of the length L is suppressed as 

C(u, A; L) ~ c L ~ D- L . (4.22) 

Therefore extended correlators vanish in the large- D limit. 
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Figure 3: The lattice in the form of a g-simplex (depicted for q = 5). 

5 The gauged Potts versus KM models 



As is known, the KM model is equivalent to the matrix model on a Bethe tree ||29|| . We 
propose in this section yet another matrix model — the gauged Potts model — which 
is equivalent to the KM model at large N providing the coordination numbers coincide. 
The gauged Potts model is convenient for studying a relation with discretized random 
surfaces and for interpreting the results of the previous sections. The proof of equivalence 
is given via loop equations which reduce at large N to a one-link equation whose different 
forms are considered. This reduction holds in the strong coupling phase where the vacuum 
expectation values of the closed Wilson loops of the gauge field vanish except for those of 
vanishing minimal area. 



5.1 The partition function 



J Potts 



(5.1) 



As is well-known P0 |, the g-state Potts model on a random lattice is equivalent to the 
matrix model 

^ x=l 

where the NxN Hermitean matrix (p x lives on the lattice which form a g-simplex (depicted 
in Fig. for q = 5). The second term in the action involves the sum over all the links (the 
link {x, y} connects the sites x and y which are nothing but the vertices of the g-simplex). 



We propose the following gauge-invariant extension of the model (|5[ 

f -it A Ntv(-y ,%)+cy < 

J UdU^fldfae V ^ ■ ^ ■ 



'GP 



(5.2) 



x>y x=l 
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where the gauge variable U xy lives on the link {x, y} and dU xy is the Haar measure on 
U(N). This construction is quite similar to the KM model [0 which is described by the 
partition function ( Q|) except the lattice is hypercubic in the latter case. 

The relation of the partition function ( |5.1Q with discretized random surfaces was stud- 
ied in Ref. |JD|. An analogous interpretation of the gauged Potts model ( |5.2| ) is based 
on the expansion (|4.2|) of the Itzykson-Zuber integral. Now the terms of the type (tr</>) 2 



generate vertices with touching surfaces |T0|. 



5.2 Loop equations 

To investigate the model (|5.2|) , let us introduce the extended open- loop averages 

G - (r "> = (s^^I^'M) (5.3) 

where T xy goes from x to y along some path on the g-simplex and the average is w.r.t. 
the same measure as in (5.2). G v \{T xy ) is symmetric in v and A due to invariance of the 
Haar measure dU under the transformation U — > U* . These quantities are quite similar 
to those (3J;) for the KM model. 

The averages ( |5.3| ) obey quantum equations of motion which are known as the 
Schwinger-Dyson or loop equations. Their derivation is analogous to the one for the 
KM model which is presented in Appendix C and the resulting equation! 



\ N v v - <j) x y X-(j)y 
9-1 / tr / 1 1 

At=1 \ IN V - <p x \~<Py 



\ N \v-<\>J N -<f> x \-<p y 

^ ( I i^f^jhr) I > • (5 - 4) 

where the sum over [i goes over 

A = g-1 (5.5) 
directions, coincides with Eq. ( |C.8| ) providing 

q-l = 2D. (5.6) 

Since A is nothing but the coordination number for the g-simplex which equals 2D for 

the hypercubic lattice, the equality Q5.6D simply means that the coordination numbers 
coincide in both cases. 



We put c = 1 in this section and in Appendix C. 
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While the loop equations which result from the variation of <fi x look similar, the cor- 
relators in the second term on the r.h.s. of Eqs. ( |5.4j ) and ( |C.8| ) are, generally speaking, 
different. This term is proportional to 8 xy which does not vanish only for x = y, i.e. only 
for a closed contour Y xx . However, the two models are equivalent in the strong coupling 
(or small-c) phase where 



\-<p x J N v A - 

= 6 oAmin (r)^(E x - K)E X + O (^) (5.7) 

in both cases as iV — > oo. The r.h.s. of Eq. (|5.7|) does not vanish only for the case of 
contractable contours T xx = with vanishing minimal area A m i n (T). Hence, the second 
term on the r.h.s. of the loop equation ( |5.4|) vanishes for T xy ^ at N = oo independently 
of whether C xy is closed or open. 

Finally, the second term on the l.h.s. of the loop equation ( |5.4| ) can be simplified at 
N = oo using the formula 

/ tr 1 1 

( n \ < t > x+iJJfi( x+ p) x ) - — - U(T xy ) - — - U Jf {r{ x+ ^) X Y xyj 

which holds providing /x does not coincide with the direction of the first link of the contour 
T xy emanating from the point x. 

The function 

OO J. oo 

F(u) = Y,F n u\ F = - (0 - £ F n r) (5.9) 

n=0 1N n=l 

which enters Eq. ([5.8|) is determined by the pair correlator of the gauge fields 



tr 

where <fi and ip play the role of external fields and t a (a = 1, . . . ,N 2 -1) stand for the 
generators of the SU(N). Eq. ( |5.10| ) is based solely on the properties of the integral over 
unitary matrices and holds [^, [31] at iV = oo. The choice of F which is not determined 
by Eq. ( |5.10| ) is a matter of convenience pj . How to calculate the function F is explained 
in the next subsection. 

Using ( |5.8| ) and introducing 

V(uo) = V'(u) - (A - l)F(u) , (5.11) 

we rewrite Eq. ( |5.4| ) and Eq. ( |C.8| ) at large N, when the factorization holds, in the same 
form 

JCi 2m [y — u>) v — X 
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where the contour C\ — the same as above — encircles counterclockwise singularities of 
the function G W \(T). 

The last term on the r.h.s. of Eq. ( |5 . 1 2| ) vanishes for T 7^ while the explicit equation 
for T = reads 

duj V(uj) _ „ ^2 



L 



Ci 2ni (A — lo) 



E u — Gi(A) = E x 



where I have denoted 



E, = 



tr 



N U 



analogously to Eq. fl2.19| ) for the KM model and defined the one-link average by 



tr 



( 



A-0 : 



■c/ f (r 



(5.13) 
(5.14) 

(5.15) 



since the r.h.s. does not depend on x and \i. Using Eq. ( |5.8| ) and noticing that 

V' (u) = V'(u) - F(u) , (5.16) 

Eq. ( |5.13| ) can be rewritten in the form of the loop equation of the Hermitean one-matrix 
model 



duj V'(oj) 



-E„, = E\ . 



(5.17) 



Id 2iri (A — to) 

The explicit equation for the case when T coincides in Eq. ( |5.12|) with one link reads 

du v'M 



G 



E v G u \ + \G V \ — E v 



(5.18) 



/Ci 2ni [y — lu) 

where G u \ is defined by Eq. (|3.7|) . Eq. (|5.13|) is nothing but the 1/A term of the expansion 
of Eq. (15381) in 1/A. 

Eq. (|5.18| ) coincides with that [|25| for the KM model providing the relation (|5.6|) holds. 
It can be shown that Eq. (|5.12| ) is satisfied for any contour V providing Eq. ( |5.18|) is sat- 



isfied. The open-loop averages ( |5.3|) for the gauged Potts model coincide with those (pT6 
for the KM model providing the lengths of the contours T xy on a g-simplex and C xy on a 
-D- dimensional lattice coincide. 



5.3 The general solution 



Eq. ( 5.18 ) looks the same as the loop equation for the Hermitean two-matrix model [13, 
[14], |j~5| . This is because at q = 2, which is associated with the Hermitean two-matrix 
model, A = 1 and the last term on the r.h.s. of Eq. ( |5.11D disappears so that one gets 
just V(u>) = V(u). To analyze it, let us consider the Hermitean two-matrix model with 
the potential 

00 a 

V($)=Y,—<p m - (5.19) 

m=l m 
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Expanding G u \ in \jv 

G« = £ + | ^ • G„(A) = ( H (fl^^Dl^)) (5.20) 
and substituting into Eq. ( |5.18| ), one gets 

Y.9mG m - i(A) = XE x -l (5.21) 

m>l 

which determines E\ versus V(A). 

The functions G n (X) are expressed via E\ using the recurrence relation 

G n+ i(\)=[ p-. 7 P^G n (u)-G {X)G n (X) , G (\)=E X (5.22) 
JCi 2,ixi (A — oj) 

which is obtained expanding Eq. (|5.18|) in 1/A. If V(A) is a polynomial of degree J, 
Eq. (|5.21|) contains E\ up to degree J and the solution is algebraic. As is proven in 
Ref. |U: 

i) Equations which appear from the next terms of the 1 / ^-expansion of Eq. (|5.18 ) are 
automatically satisfied as a consequence of Eqs. ( |5.21| ) and fl5.22| ) . 

ii) G v \ is symmetric in v and A for any solution of Eq. ( |5.21|) . The symmetry require- 
ment can be used directly to determine E\ alternatively to Eq. ( |5.21| ). 

The approach based on Eq. (|5.18|) is equivalent ^5] to that of Ref. |§ which is based on 
the Riemann-Hilbert method. To show this one takes the continuous and discontinuous 
in v parts of G v \ across the cut (cuts) which are defined by Eqs. (|3.9|) . The discontinuous 
part of Eq. ( |5.18| ) then reads 

Cont„ G v \ = 1 — — ( A + Cont^ E v - V(u)) Disc u G uX , for v G cut (5.23) 
Disc„ E v v ' 

which coincides with the equation of Ref. 0. 

To obtain a formal solution to Eq. (|5.23|) for G v \ versus E u , one notices that for any 
real v 

1 - G(z, +i o)A _ 1 - Cont„ G v \ - i Disc,, G uX _ A + Re E v - i Im E v - Viy) 



1 - G(v- i0 )\ 1 - Cont^ G v \ + i Disc„ G v \ A + Re E v + i Im E u - V'{u) 
since ~Disc u G u \ cancels at the cut (cuts) due to Eq. ( f5.23j ). 

The solution to the Riemann-Hilbert problem ( |5.24| ) for G v \ can be expressed via E\ 



as follows §, §, | 



G„ A = l-exp{-jf 



2m (u — u 



log(A-r + H) (5.25) 



where 



V(A) + F(A) , (V(\)-E X 
r ± (A) = g ±,7rp(A) = | F(A)+ ^ . (5.26) 

Eq. (|5.25| ) solves Eq. ( |5.24| ) and holds for A outside of the cut where the asymptotic 
expansion in 1/A exists. It solves, hence, the recurrence relation ( |5.22j ). 
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5.4 Alternative one-link equations 



We have obtained in the previous subsection Eq. ( [5 . 2 3| ) from Eq. ( |5.18| ). In order to show 
that Eq. ( |5.18j ) can be inferred, in turn, from Eq. ( |5.23| ), let us consider the following 
auxiliary identity 



du Cont w (E^G 



EyG 



du Cont^, G^a 



Ci 2m' v — u Jcx 2ni v — u 

du (A+ CoiA u E u -V{w)) n 



= E V G V \ — E U + / 

E V G V \ — E v + \G U \ — [ - — :-; r 

Jd 2m [y — u) 



G, 



V — U 

du Cont u {E u G wX ) 



(5.27) 

iCi 2m v — u 

which is based only on the analytic properties and Eq. (|5.23 ). Canceling the l.h.s. with 
the last term on the r.h.s., one arrives at Eq. ( |5.18| ). 

Starting from Eq. ( 5.23|) , one can derive also slightly different form of Eq. ( |5.18| ). Let 
us consider the inverse quantity 

G u \ 



It is convenient to introduce 
and 

so that 



G u x - r , 

1 _ ^uX 

%x = Gvx — 1 
%x = G u \ + 1 



1 — G u \ 



Tux — %\ ■ 
Eq. ( |5.18| ) can be written in terms of T v \ as 

du V'(uj 



Ci 2m {y — uj) 



%jX — E v %x + XT uX + A , 



(5.28) 

(5.29) 
(5.30) 

(5.31) 
(5.32) 



It follows from Eq. (|5.31 ) that 

Disc^ T uX 
Cont^ f uX 

so that Eq. ( p.23| ) can be rewritten as 

1 



Disc,, 77, 



v J-vX 



\%x\ 2 ' 
Cont^ %x 

\ vX\ 



(5.33) 



Cont,, 77, 



v 2 uX 



Disc„ E v 



A - Cont,, E v - F(u) Disc,, % x , for v G cut . (5.34) 



The solution to the Riemann-Hilbert problem (|5.34|) which is analogous to ( [5.25| ) reads 

duj 1 



71/ a = exp 



Ci 2m {y — uj) 



log (A - r-(u)) 



(5.35) 
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Moreover, these two solution coincide due to the relation ( |5.31|) . 

The condition for the r.h.s. of Eq. ( |5.25|) (or (|5.35|) ) to be symmetric in v and A is 
given H by Eq. ( |3.14| ). Since G v \ is symmetric in v and A, the master field equation || 

ICi 



(5.36) 



obtained as the \jv term of Eq. fl5.25| ), will be satisfied as a consequence of Eq. ( |3.14j ) 
which guarantees the symmetry. 

To derive an analogue of Eq. (|5.18|) starting from Eq. ( |5.36| ), one proceeds similarly 
to ( CT : 



did Cont w (E U T U , 



c x 2m 



V — Id 



E V T U> 



did E„, Cont,,, % 



uX 



c x 2m 



V — Id 



EJT V \ + 



did 



E v %\ 
Hence, we get 

which yields 



\% x + A + 



Cx 2m 
du F{lo) ~ 

*uX + 



-A+ Cont^ E^ + F {id)) ~ 



V — id 

duo Cont,, EjT^x 



Cx 2m {y — id) 

did F(id 



Cx 2m 



v — Id 



Cx 2m [y — id 

did F(id 



-E v T v \ + XT U \ — A 



G, 



LdX 



—E u G v \ + AG 



vX 



(5.37) 



(5.38) 



(5.39) 



iCx 2m [y — id) 

which is an analog of Eq. ( |5.18| ) with V' replaced by F. Note that with the relation ( |5.28| ) 
Eq. ( |5.39| ) differs from Eq. ( p,18| ) only by the sign of the first term on the r.h.s.. 

Eqs. Q5.21 ) and ( 5.40 ) are equivalent. It is a matter of practical convenience which 
equation to solve. Eq. ( |5.18| ) is more convenient when V'(id) is a polynomial when it 
reduces to an algebraic equation for E v . Eq. (|5.39| ) in turn is more convenient if F(u) is 
a polynomial when it reduces to an algebraic equation for E v 



F m G m -i(\) — \E\ — 1 

m>l 

and G m (\), which are defined via the asymptotic expansion 

G v \ 



(5.40) 



V 



n=l 



V 



n+1 



(5.41) 



can be expressed in terms of E\ using the recurrence relation 

did F(id] 



G n +i(\) 



Cx 2m (A — to) 



G n (u) + G (\) G n (X) 



G (X) = E X . 



(5.42) 
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The explicit solutions exist at any A for the quadratic potential || and the logarithmic 
potential ( |2.1| ) when one gets [7] 

a + 1 



V(A) 



a 



F(X) 



(5.43) 



b — A a + A 

It is easy to verify that Eq. ( |5.39| ) is satisfied for G u \ given by ( |5.28| ) and ( |3.10| ) analogously 
to Eq. (KW). 



6 Discussion 

The KM model with the logarithmic potential ( j2.1| ) is a very nice, explicitly solvable 
example of how a critical behavior can exist for D > 1. The analytic properties of the 
solution and the position of the eigenvalue support are for a < of the type advocated in 
Ref. 0. From this point of view it illustrates how such unusual analytic properties can 
emerge in matrix models. 

The continuum theories which are obtained in the vicinity of the critical lines are asso- 
ciated with the jstr = and •y str = —1/2 phases of 2D gravity with matter independently 
of the fact that one starts from the matrix model on a D- dimensional lattice. Similar 
results can be obtained for D — 1 if one considers a closed matrix chain (for an open one 
the gauge field can be absorbed by a gauge transformation of <ft x ). Therefore, an exact 
solution to this D = 1 problem is found for the potential ( |2 . 1| ) . 

While the phase with 7 str = seems to coincide with the standard d = 1 string, a 
Kosterlitz-Thouless phase transition which occurs at the tri-critical point separates it 
from the phase with 7 str = —1/2 which is, however, a novel one since the continuum 
limit of matter at large distances sets up due to a special behavior of the Itzykson-Zuber 
correlator of the gauge fields. This phase never shows up, say, in the case of an open 
matrix chain where the gauge field can be gauged away. 

The existence of these two continuum limits does not mean that a d > 1 phase is 
impossible for the KM model with the logarithmic potential on a D- dimensional lattice. 
An example is a two-matrix model where a nontrivial phase with 7 sir = — 1/3 realizes at a 
tri-critical point while generically it has 7 str = —1/2. It is most interesting to investigate 
the KM model with the potential (|2.1|) in the vicinity of the tri-critical point which is not 
done in the present paper. The point is that a singular behavior of the Itzykson-Zuber 



correlator, which occurs at the tri-critical point, may induce [32] a phase transition to a 
d > 1 stringy phase. 

It could be, however, that the logarithmic potential ( |2.1| ) is too simple to exhibit a 
nontrivial scaling behavior except for 7 str = and 7 s t r = —1/2. From this point of view 
the results of Sect. 4 about 7 s t r > in the large- .D limit for other potentials, including 
the quartic one, may be interesting. I do not think, however, that an exact solution can 
be obtained for these potentials. One should look rather for a scaling solution just near 
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the critical point. In particular, l/Z)-corrections is the simplest problem to study. I hope 



that the fact that the mechanism of Ref. [IIJ of getting ^ str > realizes for the KM model 
in the large- .D limit can answer the question whether it corresponds to branch polymer 
or stringy phases. 
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Appendix A The eigenvalue support for a < 

To find the position of the contour C in the complex A plane along which the spectral 



density, p(A), has the support for a < 0, we apply the criterion f2C§ which is based for 
the one-cut solution (p. 22 ) on the following integral 



G(X) 



V\t) - 2E T 



(A.l) 



The eigenvalue support is along the curve from x _ to x + determined by 

ReG(A)=0 for A G C (A.2) 

which is embedded into the domain where ReG(A) < 0. The spectral density is given by 

1 dG(\) 



P(A) 



for AeC. 



2wi d\ 

Let us consider in some detail the case a = — 1 when 

z — f3 , X— — X-]- — f3 — 



1 



(A.3) 



(A.4) 



according to Eqs. ( |2.56| ) and ( |2.37| ). Since the points x_ and x + coincide, C forms a closed 
loop which passes through the point (3 — 1/213 and encircles the point (3. 

The one-cut solution ( |2.33| ) has for a = — 1 the following explicit form 



2(X-(3) 



2(X-(3) 




(A.5) 



where the plus sign should be substituted when A is outside of the loop C and the minus 
sign is associated with A which is inside C. The expression (|A.5| ) has the same expansion 
in 1/A as an analytic continuation of the solution for positive a to a = — 1 which were 
give for E\ simply 1/(A — (3) with the pole at A = j3. The true solution ( |A.5|) is analytic 
at A = (3. 
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Figure 4: The eigenvalue support for a = — 1 and (3 = 1/2: the loop to the right of 
Re A = —1/2. The whole curve represents G(X) = 0. 



Substituting (|A.5|) in to Eq. (|A.1| ) one gets 



G(A) = ± -2/3A - log (A - 0) + 2f3 2 - 1 - log 2/3 



Introducing the new variable k by 



(A.6) 



(A.7) 



1 , ( 1 



which coincides with A at /3 = 1/2, one rewrites Eq. ( |A.6|) in the form 

G(k) = T 

which does not depend explicitly on (3. Denoting 

Re k = u , Im k = v , 
we find the solution to the equation Re G(k) = to be 



-2m- 



(A.8) 



(A.9) 



(A.10) 



The curve (|A.10| ) is depicted in Fig. ^ and is quite similar to the one for the Penner 
model ||18|| . The contour C which represents the eigenvalue support is the closed loop to 
the right or u = — 1/2. The lines to the left of u = — 1/2 separate regions of positive and 
negative G(X) which is positive to the left and negative to the right of these lines. Since 
k = A at j3 = 1/2, Fig. | represents the eigenvalue support for the (3 = 1/2 while that for 
an arbitrary (3 > 1/2 can be restored using Eq. ( |A.7|) . 
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The variables (|A.9[) are convenient to perform calculations with the spectral density 
which is determined by Eqs. ( |A.5| ) and ( |A.10[ ) to be 



p(A) = ^ f e 2u+ \u - ± - iv) + l} . (A.ll) 
While this expression is complex, the element of probability 



e -"- 1 / 2 4- 2v 

dXp(X) = du . — (A. 12) 



with v given by Eq. ( |A.1(J| ) is real and positive. It is easy to verify the normalization 
condition 

"* , e~ u - 1/2 + 2u , k . 

du — j = 1 (A.13) 

1/2 Je-^-(u 1X2 



2 

where = 0.77846454 is the solution of the equation 



e-«.-V2_ u#+ 1 =0 . (A.14) 



Appendix B The convolution formula 



We prove in this appendix that ( |3.22| ) satisfies Eq. (|3.23| ) 



Let us first change the variables x — >l/x, y — >l/y and rewrite Eq. ( |3.22j ) as 

C c (-,-;Vu~)=x 2 y 2 C(x,y;Vu~) (B.l) 
x y 

with 

d ( X >V>^ = ( X -y)2 + Mx + y)+U*- (R2) 

Substituting t — > 1/t, we rewrite Eq. ( |3.23|) as 

- dty/iC(x,t;y/u)C(t,y;y/v) = C(x,y; ^1 + V^) ■ (B.3) 

7T JO 

To prove Eq. ( B.3|) , let us calculate the integral on the l.h.s. taking residues at t = t±(x) 
and t = t±(y) with 

t±(x) = x — u± 2i\[xu (B.4) 

and 

t±(y)=y-v±2i^yv (B.5) 

being the poles of C(x, t; ^Ju) or C(t, y\ y/v), respectively. Since t+(x) and t-(x) (or 

and t-(y)) lie on the opposite sides of the branch cut of \fi which goes along positive real 

axis as in depicted in Fig. j|, we get 
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© 



t+(x) t+(y) 

• • oo 



Figure 5: The position of the cut of p c {t) in the t-plane (the bold line) and the positions 
of the poles of the integrand in Eq. (|B.9| ), t±(x) and t±(y), which are given by 
Eqs. flB.4|) and ([B.5D, respectively. 



\Jt + (x) = yfx + iy/u , \Jt-{x) = —y/x + iyfu (B.6) 



and 



yjt+(y) = y/y + iy/v , yjt_(y) = -y/y + lyfr . (B.7) 
Collecting the contribution from all four poles, factoring the common denominator as 

[(x — y) 2 + 2(x + y)(u + v) + (u — v) 2 } 2 — 64xyuv 
= [(x - yf + 2(^5 + v^) 2 (x + y) + (V^ + v^)1 

■ [(x - y) 2 + 2(V5 - v^) 2 (^ + y) + (>/u - v^)1 (B.8) 
and cancelling the last factor with the numerator, we get finally 
1 dty/i 



it Jo [(x - t) 2 + 2u(x + t) + u 2 } [(y - t) 2 + 2v(y + t) + v 2 } 
fu + y/v 1 



2^ \{x-y) 2 + 2{y/^+^) 2 {x + y) + {y/^+y^YY (B ' 9) 
This completes the proof of Eq. QB.3|) . 

Appendix C Derivation of the loop equations 

Loop equations of the KM model or the gauged Potts model result, as usual, from the in- 
variance of the integration measure under an infinitesimal shift of fields. One distinguishes 
equations resulting from the shift of the matter fields <p x (that gives, in particular, the 
lattice Klein-Gordon equation) and that of the gauge field U xy (that gives, in particular, 
the lattice Maxwell equation in the case of the standard lattice gauge theory). 
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C.l Matter loop equation 

Let us consider an equation which results from the invariance of the measure over in 
the open-loop average ( |3.6| ) under an infinitesimal shift 

+ (CI) 

of 4> x at the given site x with £ x being an infinitesimal Hermitean matrix. For <f> x being a 
matrix from the adjoint representation of SU(N), one should impose tr^ x = in order 
for the shifted matrix to belong to the adjoint representation as well. For the general 
Hermitean matrices <f> x , £ x is arbitrary Hermitean. 

These two cases can be considered simultaneously introducing N 2 generators 

[t% = [t%) (C.2) 

with t a (a = 1, . . . , N 2 -l) being the standard generators of SU(N). The generators (|C.2|) 
obey the following normalization 

t It A t B = 5 AB (C.3) 

and completeness condition 

[t A ] ij [t A ]ki = NS i iS kj . (C.4) 
An arbitrary N x N Hermitean matrix can be represented as 

= t V where <j> A = ( ^<f>, | t a cf) (C.5) 



with = vanishing if is taken in the adjoint representation of SU(N). 

To derive the loop equation I apply a trick similar to that used in deriving loop 
equations of QCD |l7j]. Let us consider the loop average 



where the averaging is taken with the same measure as in Eq. fll.lj ), which vanishes due 
to the gauge invariance. Performing the shift ( CI ) of <f> x , using the invariance of the 
measure and calculating d/d<f) x , one gets 



D 



= / ^(t A ^—t B ^—u{c xy )-^—u\c xv 
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91 



X 



x + /J, 



a) 



b) 



Figure 6: The graphic representation for G v \(C xy ) (a) and G v \(Cr x+l A x C xy ) (b) entering 
Eq. ( |5.12 ). The empty circle represents 4> x +ii while the filled ones represent 
1 The oriented solid lines represent the path-ordered products 
The color indices are contracted according to the 



or 



f — qJx \~4>y ' 

U{C xy ) and U(C( x+ ^ x C xy ) 
arrows. 



The l.h.s. of this equation results from the variation of the action while the r.h.s. represents 
the commutator term resulting from the variation of the integrand. 

The averaging over the gauge group picks up two nonvanising invariant equations for 
the Hermitean matrices. The first one can be obtained contracting Eq. ( U77) by 5 AB 
(A, B = 0, . . . , N 2 -l) while the second one is given by the A, B = component. 

The first equation for the Hermitean model reads 




U(C( x+ll ) x 





tr 



1 



tr 



1 



NV- 



-U(G 



xy, 



A - 



+ S X y 



tr 



N \u 



(^—U(C xy )-^—) - f-J_£/t(c 




(C.8) 



where the path C( x+fl ) x C xy on the l.h.s. is obtained by attaching the link (x, x + fi) to 
the path C xy at the end point x as is depicted in Fig. || Using the definition ( p.6|) , this 
equation can be written finally at large iV in the form ( |5.12|) . 
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The second equation which is given by the A, B = component of Eq. ( |C.7|) reads 
tr . x A tr \ tr 1 




(C.9) 

In the large- N limit when the factorization holds, Eq. ( |C.9| ) is automatically satisfied as 
a consequence of the 0((VA) _1 ) term in Eq. 



When the matrices <fi x belong to the adjoint representation of SU(N), there exists 
only one invariant which results from contracting Eq. (|C.7|) by 5 ab (a,b = 1, . . . , N 2 -l). 
The point is that the B = component does not appear since solely the variation d/d(p b x 
is permitted for the adjoint-representation matrices. Therefore, in order to obtain an 
analog of Eq. ( |C.8| ) for the adjoint scalars, one should subtract the l.h.s. of Eq. (|C.9| ) 
from its l.h.s. and the r.h.s. of Eq. ( |C9[ ) from the r.h.s. to kill the A=B=0 component. 
The result differs from Eq. ((J. 8) by contact terms which enter Eq. ( |U.9| ) and should not 



survive as N — > oo. The vanishing of these contact terms in the large- N limit, when the 
factorization holds, is obvious for an even potential V{4>) = V{— (/)). 



C.2 Maxwell open-loop equation 



The open-loop averages ( |3.6| ) obey one more loop equation which results from the imparl- 
ance of the Haar measure over U x r x+t A under the shift 

U x {x+n) — ► (1 + itx(x+ii))Ux{x+ij) (C.10) 
of U x ( x+fl ) at the link (x, fi) with e x ( x+tl ) being an infinitesimal traceless Hermitean matrix. 



The loop equation can again be obtained by the trick p7|| . Let us consider the loop 
average 



= 0, (C.ll) 

where the averaging is taken with the same measure as in Eq. Ql.l] ), which vanishes due 
to the gauge invariance. Performing the shift flC.10|) of U z [ z+I1 ) at some link (z, ji) G C xy 
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t c z 
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Figure 7: The contours C xz (C zx ) or C zy (C yz ) (a) and C^^ x {C x{x+ ^) (b) for 
Eq. ( pl2|) . 



and using the invariance of the Haar measure, one gets 



tr 

N V 



N 



(C.12) 



where the contours are depicted in Fig. |7|. 

The desired loop equation can be obtained contracting Eq. ( |C.12| ) by the structure 
constant f abc which is defined by the commutator 



[t a i t b ] =if abc t c . 
The normalizations are fixed by Eq. ( |C.3| ) which gives 

jabc jdbc 2,N5 ab 

I shall utilize as well the following formulas 

if abc [t a }ij[t b }ki = [t c }uhj - [t c ]kj$u 

and 

ir c [t a ] ij [t%[t c ] mn = n(s u s 

kn^mj ^in^kj^ml) 



(C.13) 



(C.14) 



(C.15) 



(C.16) 



Eq. ( |U.16|) can be obtained from Eq. (|C.16|) multiplying by [t c ] mn and substituting the 
completeness condition Eq. ( |C.4|) . Eq. ( |C. 16|) was used, in particular, in studies of 
loop equations of QCD. 
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«z 



Figure 8: The graphic representation of the r.h.s. of Eq. ( |C • 17| ) • 



The result of contracting Eq. ( |C.12| ) by %f ahc reads 



tr ' ' U{C xy )-±-U\C zy 



\ N \u- 
1 / tr / 1 



tr 1 tr 1 



2 \ N V- 



A - <f)y 

U (C XZ ))4> Z U (C z ( z+Il ))(f) z+ ^U (C( z+fl )y 



N v-(f) x N \ -<f) y / 



-U(C, {l+r) )^ z+ ^{C, {I+ ^,U(C, v )- -U\C xy ) 

A — <py 



+u(c xy ) 
-u(c xy )- 



A - 4>y 

1 



U ] (C [ z+ljL )y)(t) z+IJ U ] {C z i yZ+li ))(t) z U ] (C x ( z+ ^ ) 



U\C zy )(j) z U(C z{z+ ^)<p z+ ^U\C x i yZ+ll )) 



(C.17) 



A - (f>y 

where the contours entering the r.h.s. are depicted in Fig. ^| 

Note that quantities of a new type appears on the r.h.s. of Eq. ( |C.17| ) so that the set 
of equations is not closed. The meaning of this equation is that it expresses the one-link 
correlator of UUU^U^ via the correlator of UU', i.e. via C(x,y). 
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